#include "leapfrogM.h"
#include "MatrixV.h"
using namespace std;
int main()
{
    double h = 0.05, T = 17;
    double k = 0.8*h;
    int m = 1/h - 1;
    vector<double> U0;
    for (int i = 0; i < 9*(m+1)+m; i++)
    {
	U0.push_back(u1((i+1-2)*h));
    }
    cout << "==============================" << endl;
    print(U0);
    cout << "==============================" << endl;
    leapfrog_method M1(h, T, k);
    cout << "==============================" << endl;
    print(M1.Solver());
    cout << endl;
    k = h;
    leapfrog_method M2(h, T, k);
    cout << "==============================" << endl;
    print(M2.Solver());
    cout << endl;
}
